###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)
require(ExpDes)


#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
source("funcoes.txt")



#########################################################################
############################Análise do DBC###############################
#########################################################################

#Mudar diretório
setwd("C:/Users/user/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2016/dados")

##Importar dados
dados=readXL("jumento.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)


Bloco=as.factor(Jumento)

######Análise de Variância#
model=aov(VR~Bloco+Diluente)
anova(model)
cv.model(model)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)



#Teste de Shapiro-Wilk - NORMALIDADE
shapiro.test(e) 
#Teste de Oneill-Mathews - Homogeneidade de Variância
VR=VR
Bloco=Bloco
Trat=Diluente
oneillmathewsDBC(VR,Diluente,Bloco)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model)

###Teste de aditividade dos efeitos de tratamento e bloco
df<-df.residual(model)
MSerror<-deviance(model)/df
nonadditivity(VR, Bloco, Trat, df, MSerror)

####Anova com Teste de Tukey####################
rbd(Diluente,Bloco, VR, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05)

####Anova com Teste de Scott-Knott####################
rbd(Diluente,Bloco, VR, quali = TRUE, mcomp = "sk", sigT = 0.05, sigF = 0.05)